The universality class of fluctuating pulled fronts 
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It has recently been proposed that fluctuating "pulled" fronts propagating into an unstable state 
should not be in the standard KPZ universality class for rough interface growth. We introduce an 
effective field equation for this class of problems, and show on the basis of it that noisy pulled fronts 
in d+1 bulk dimensions should be in the universality class of the (d+lJ + lD KPZ equation rather 
than of the d+TD KPZ equation. Our scenario ties together a number of heretofore unexplained 
observations in the literature, and is supported by previous numerical results. 

PACS numbers: 5.40+j, 5.70.Ln, 61.50.Cj 

Consider spatio-temporal systems in which the impor- 
tant dynamics is governed by the propagation of fronts 
or interfacial zones separating two domains whose bulk 
dynamics is relatively trivial or uninteresting. In the 
presence of fluctuations, the theory of the stochastic be- 
havior of such fronts or interfaces is well-developed . 
In particular it is known that many such fluctuating d- 
dimensional interfaces in d+1 bulk dimensions are de- 
scribed by the KPZ equation jsj for their height h, 
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with rj a random gaussian noise with correlations 
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We will follow common practice to refer to this equa- 
tion as the d-f _ZD(imensional) KPZ equation, where the 
d refers to the dimension of the interface and the +1 io 
the time dimension; denotes the coordinates perpen- 
dicular to the direction of propagation of the interface. 

The fact that the scaling behavior of so many stochas- 
tic interfaces fall in the d+lD KPZ universality class, is 
due to the fact that (1) this equation contains all the 
terms in a gradient expansion which are relevant in a RG 
sense; and (2) that the long wavelength deterministic dy- 
namics of many interfaces is local in space and time, i.e., 
of the form u„ = w„(V/i, V^/i, • • •), expressing that the 
normal velocity Vn becomes essentially a function of the 
instantaneous slope (angle) and curvature of the interface 
only. Upon expanding in the gradients, adding noise, and 
retaining only relevant terms, one then arrives at (Q). 

The starting point of such an argument, the fact that 
one can integrate out the internal structure of the inter- 
face and on long length and time scales think of it as a 
mathematically sharp boundary with effective dynamics 
expressed by a boundary condition w„=u„(V/i, V^/i, • • •) 
which is local in space and time, is appealing and usually 



correct. Intuitively, one associates it with the interfacial 
zone being sufficiently sharp on a spatial scale. Nev- 
ertheless, there have been scattered observations in the 
literature which indicate that there is more to it: (a) 
Some continuum reaction-diffusion equations have prop- 
agating planar interfaces of finite width which are stable, 
but which become weakly unstable for discrete particle 
model equivalents |4| , contrary to what the above coarse- 
graining picture would suggest, (b) The empirical rela- 
tion observed for the distribution of DLA fingers in a 
channel and the interface shape of a viscous finger could 
not be understood from the standard continuum model 
until the innocuously looking reaction term was regu- 
larized [^; on hindsight, this was because the standard 
mean-field DLA equations do not give the appropriate 
"local" boundary conditions of the type w„ — — /iV„p. 
(c) In a simple stochastic particle model with fluctuat- 
ing fronts, non-KPZ scaling was observed contrary to 
what one would naively have expected. 

It turns out that these observations all have one com- 
mon denominator [@,||, in that they are related to the 
existence of two classes of fronts, "pushed" and "pulled" 
fronts. Pushed fronts are the usual ones: their dynam- 
ics is determined by the behavior in the interfacial zone, 
a region of finite thickness, and their response to the 
bulk fields is local in space and time Pulled 
fronts, on the other hand, propagate into a linearly un- 
stable state. Although they do not differ noticeably 
from pushed fronts in their appearance, their dynamics 
is driven by the growth and spreading of perturbations 
about the unstable state in the semi-infinite region ahead 
of the front [^; hence they are particularly sensitive to 
slight changes in the dynamics there |TT|J^ . These impor- 
tant differences led two of us |^ to propose recently that 
fiuctuating variants of d-dimensional pulled fronts in d+1 
bulk dimensions would indeed not be in the d+lD KPZ 
universality class, even though pushed fronts do effec- 
tively give local boundary conditions on long length and 
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time scales, and hence do give rise to d+lD KPZ scaling 
in the absence of coupling to a diffusion or laplace field in 
the bulk . Simulations of a simple stochastic lattice 
model were consistent with these arguments, and with 
the earlier observations of 

In this paper, we will argue that fluctuating pulled 
fronts are indeed in a different universality class from 
the usual pushed ones which show the standard KPZ 
behavior. Indeed, we will show that the semi-infinite re- 
gion ahead of the front can not be integrated out, and 
effectively enhances the dimension by i: we introduce 
a field equation for fluctuating pulled fronts and argue 
that d-dimensional fronts in d+1 bulk dimension are in 
the universality class of the (d+l)+lD KPZ rather than 
the d+lT) KPZ equation. This surprising scenario, which 
also builds on the insight of for the stochastic behav- 
ior of pulled fronts in one bulk dimension, is fully con- 
sistent with our earlier 2D simulations and also with 
the heretofore unexplained results of in higher dimen- 
sions. In addition, as we shall discuss, our scenario leads 
to a number of interesting new questions and challenges. 

A stochastic equation for pulled fronts should obey 
two requirements: in the usual stochastic lattice mod- 
els with fronts, no particles are spontaneously generated 
when there are none already. Secondly, the average front 
speed and the local fluctuations ahead of the front re- 
main always flnite in such lattice models The field 
equation should be consistent with these basic facts. So 
when we consider a stochastic field equation for 4'{x, rj_, t 
in d+1 dimensions (x,r^) of the type 
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these requirements put constraints on the function / and 
the noise term g{(t>)'ri. The stochastic noise ri{x,r±,t) 
has delta correlations as in (^), and is interpreted in the 
Stratonovich sense, but our arguments will not rely on 
the distinction between Ito and Stratonovich calculus. 
Stochastic field equations of this type have e.g. been 
used already since long for studying the scaling behav- 
ior of homogeneous bulk phases like directed percolation 
1^^ ; investigations of noisy fronts in such equations are 
more recent — see e.g. for an analysis of stochas- 
tic pulled fronts and a discussion of the applicability to 
various systems. Here we focus on the proper form for 
an effective stochastic field equation for pulled fronts. 
For /(0), which determines the dynamics of determinis- 
tic fronts in the absence of noise, we choose the standard 
form for pulled front propagation f{(j>) — (p — (j)^, which 
gives saturation of the field (j) behind the front where 
0^1. How should the noise term g{(j>)T] look like p7|? 
The requirement that if there are no particles {(t>=0) none 
are created spontaneously, implies that there should be 
no additive noise term, and hence that g((j)—0) — 0. For 
(j) nonzero but small, it is natural to assume a power 
law behavior g{(j)) ~ (p"; in the studies of the homoge- 



neous bulk properties of directed percolation the choice 
a = 1/2 was been made [^Sj, motivated by the idea that 
typical bulk fiuctuations are of the order of the square 
root of the particle density. For pulled fronts, however, 
the dynamically important region is ahead of the front, 
where </) — > 0. Our second requirement that the rela- 
tive fluctuations g{(j))rj/4' remain finite here shows that 
the natural choice is a = 1, i.e., g{4') (j) for <C 1 
p8[ . The linearity of g for small (j) is sufficient for our 
subsequent analysis. In our numerical studies, we have 
actually taken g{(j>) = (j>{l — (p) , a form taken to suppress 
fluctuations behind the front. This makes it numerically 
easier to focus on the fluctuations of the front position 
itself, without affecting the essential results. Specifically, 
we thus propose as the generic stochastic field equation 
for pulled fronts 
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Let us now turn to the analysis of the stochastic be- 
havior of fronts which propagate along the x-direction 
into the linearly unstable state 0=0. The crucial feature 
of pulled fronts is that even though the full dynamics of 
the fronts is nonlinear, it is essentially determined in the 
"leading edge" , the region ahead of the front where re- 
mains small enough that the nonlinear saturation term 
—(j)^ which limits the growth, plays no role: the linear 
spreading and growth of perturbations about the state 
0=0 almost literally "pull the front along". An impor- 
tant recent development has been the realization that 
this simple intuitive picture can be turned into a sys- 
tematic scheme to calculate even the convergence of the 
front speed to its asymptotic value v* . Remarkably, this 
relaxation is governed by universal power laws which can 
be calculated exactly even for general equations ^ . The 
fact that the stochastic fiuctuation effects that we want to 
investigate are dominant relative to the deterministic ve- 
locity relaxation terms, suggests to calculate these along 
similar lines. For the deterministic case (77=0), fronts in 
(^) propagate with an asymptotic speed v* = 2y/D. In 
a frame ^ = a; — v*t moving with this speed, the asymp- 
totic front solution has an exponential fall-off ^ 



with A* = 1/v-D for large positive ^. The asymptotic 
relaxation analysis of deterministic fronts is based on the 
so-called leading edge transformation = ^ipi which 
transforms into 



^ = DV'^^p + - (1 + ri)^^e' 



(6) 



Here ip has been written in the frame moving with veloc- 
ity V* in the x-direction: ip — , r_L, t). 

In the analysis of deterministic fronts |^ , the nonlinear 
term on the right hand side (which is exponentially small 
for ^ ^ 1) essentially plays the role of a boundary condi- 
tion for the semi-infinite leading edge region where is 
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small — it allows the nonlinear region to match properly 
to the leading edge which "pulls" the front. As explained 
above, this holds a fortiori for fluctuating pulled fronts: 
their stochastic fluctuations are essentially determined 
by the region where the linearized equation can be used. 
Now, as is well known, upon making a Cole-Hopf trans- 
formation = e^, the linearized equation transforms to 



dh 
'dt 



= DV^h + D{Vhf 



V 
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which is nothing but the (d+l)+lT) KPZ equation 
for the d+1 dimensional field /i(^,r^,i)! 




FIG. 1. Left panel: snapshot of the field 4> at time i = 20 in 
a 2D simulation of (j^ with D = 1 and e — 10. The thick line 
is the position of the front, defined by tracking the line where 
(j){£,,r^,t) = 1/2. Right panel: the same data as in the left 
figure, plotted in terms of the height variable h. Note that h 
has the appearance of a (slanted) fluctuating surface. The flat 
portion on the left is the region behind the front and where 
h ~ A*^ since (f> ^ 1. The thick line indicates the height fluc- 
tuations along a line of constant ^. This illustrates that the 
one- dimensional position fluctuations along the pulled front 
illustrated by the thick line in the left panel are related to the 
height fluctuations of the two-dimensional fluctuating surface 
of the leading edge variable h. The scaling behavior of these 
is that of the 2+lD KPZ universality class. 

As illustrated for two bulk dimensions in Fig. |l|, the ID 
fluctuations in the front position in the propagation di- 
rection arc defined by tracing a line where (j) ~ const ^ e.g., 
(f) = 1/2. Since (p ~ e^^ the front fluctuations in the 

direction are given by £,{r±,t) = h{£^,r±,t)/ X* -I- ~ 
h{£_o, r±,t)/\* + £_o, where the constant is determined 
by the level curve of which we trace to determine the 
front position. Thus indeed the position fluctuations of a 
d-dimensional pulled front in d+1 bulk dimensions map 
onto the height fluctuations along a line of a KPZ surface 
in d+1 dimensions — see Fig. ^ The growth and rough- 
ness exponents are therefore those of the (d+l)+lT> KPZ 
universality class! 

The above scenario unifies a number of different re- 
sults. It can immediately be compared with the simu- 
lation results of the stochastic lattice model of In 
that paper a 2D lattice model was introduced in which 
by changing a simple birth and death rule of particles 
ID fronts could be tuned from pushed to pulled. The 
scaling exponents of the pushed model were found to be 
the standard 1+lD KPZ ones, as it should, while those 



of the pulled variants were close to those of the 2+ ID 
KPZ universality class. More importantly, without any 
adjustable parameters, the distribution functions for the 
long-time saturated width of the fronts in this model for 
finite transverse width L±^ |19| are completely in accord 
with our scenario Moreover, although fronts in ID 
do not have transverse fiuctuations, the wandering of the 
position of pulled fronts in one dimenension is also con- 
sistent with 1+lD KPZ scaling |l^. Finally, the observa- 
tions of Riordan et al. [|| that in three (and higher) bulk 
dimensions their fronts did not appear to show a power 
law growth of the front width finds a natural explanation. 
According to j|] their fronts are pulled and so they should 
be governed by the 3+lD KPZ equation. The free A = 
fixed point in this equation is stable and has no diver- 
gent interface width. Apparently above two dimensions 
the model of |^,|| is in the weak-coupling limit. 

On hindsight, our arguments also justify the regular- 
ization of of the mean-field equations for DLA in a 
channel: the full problem involves pushed fronts but the 
mean-field equations have pulled front solutions. The 
regularization effectively cures this by making the fronts 
into pushed ones. 




FIG. 2. The increase of the root mean square front width 
W = [([/i(rx,t) - 7i(rx,t)]2)]^/2 (with the overbar denoting 
an average over rx) as a function of time. Data are for simula- 
tions of Eq. (^ both with nonlinearity (full line) and without 
nonlinearity (dashed line), for e = 5 and an effective diffusion 
constant D = 0.4, which corresponds to dimensionless KPZ 
coupling constant A = 25 The front position h in the x 
direction is defined as the level line where (j> = 0.5. The fact 
that the growth exponent is essentially the same with and 
without nonfinearity in the (/> equation justifies our assertion 
that these terms do not affect the dominant scaling behavior 
of pulled fronts. 

The validity of the crucial step of our derivation, the 
assertion that the nonlinearities in (j|) or (^) can be ne- 
glected because the leading edge where is the es- 
sential region, can be tested independently. In Fig. ^ we 
show simulation data of the wandering of the lines where 
(f> ~ 0.5 in Eq. (^ in 2D, both with and without the 
nonlinearity. Following |E0|, where the linearized version 
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of was already employed to study the 2+lT) KPZ ex- 
ponents numerically, we have taken parameters so as to 
make the dimensionless coupling A = D^tjv^ k, 25. This 
value appears to be close to the fixed point value and so 
slow transients are minimized |20|. We find indeed that 
the two datasets with and without the nonlinear term in 
show the same growth exponent, with a value close 
to the one /3 « 0.24 of the 2+lT) KPZ equation. This 
gives confidence in the validity of our assertion that the 
nonlinear terms in the front equation are not important 
for the scaling behavior of pulled fronts. 

The main steps of our line of argument are elegantly 
direct and build on various previously established ideas; 
at the same time our scenario also raises a number of new 
questions and challenges for further research: 
( i) There is no systematic theory for the transition from 
the pushed to the pulled regime in stochastic lattice mod- 
els, so it is difficult to determine a priori which models 
lead to the standard pushed fronts and which ones to 
pulled ones. E.g., fronts in the directed percolation prob- 
lem are pushed and obey KPZ scaling in one special case 
p2t , but it is not known whether this is generally so. 
(a) Finite size scaling of the KPZ equation is nor- 
mally done for interfaces of size L± in all directions. 
Our scenario, on the other hand, leads one to con- 
sider anisotropic scaling, since there is effectively a time- 
dependent cutoff in the ^-direction ||l3|. The crossover 
scaling is completely unexplored, but is most likely quite 
tricky: for fixed L± the results of |13| for fronts in one 



dimension suggest that one should see subdiffusive wan- 
dering of the average front position, ((/i)^) ~ Vi (rather 
than ^ t) because the cutoff in the ^-direction grows as 
^ \fi, but our simulations seem to suggest that the 
crossover to this regime happens at such extremely long 
times that it can not convincingly be seen in practice. 
Moreover, the crossover is likely to depend significantly 
on the initial conditions fl^ . 

(Hi) According to the results of [pP,|7|], pulled fronts are 
very sensitive to finite particle effects, so that the con- 
vergence to a continuum limit is extremely slow. This 
crossover is poorly understood, especially in higher di- 
mensions where it is not even known how important the 
effects are in practice. 

In conclusion, we have put forward an effective field 
equation for pulled fronts and argued on the basis of 
it that pulled fronts in d+1 bulk dimensions are in the 
(d+l)+lD KPZ universality class rather than the d+ID 
KPZ universality class. The scenario ties together var- 
ious results in the literature and brings up various new 
issues for future research. 
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